Isospin dependence of 6 He+p optical potential and the symmetry energy 
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A consistent folding analysis of the elastic p( 6 He, 6 He)p scattering and charge exchange 
p( 6 He, 6 Li*)n reaction data measured at Ei^b = 41. 6A MeV has been performed within the coupled 
channels formalism. We have used the isovector coupling to link the isospin dependence of 6 He+p 
optical potential to the cross section of p( 6 He, 6 Li*)n reaction exciting the + isobaric analog state 
(IAS) at 3.563 MeV in 6 Li. Based on these results and the Hartree-Fock calculation of asymmet- 
ric nuclear matter using the same isospin-dependent effective nucleon-nucleon interaction, we were 
able to confirm that the most realistic value of the symmetry energy E sym is around 31 MeV. Our 
analysis has also shown that the measured charge exchange p( 6 He, 6 Li*)n data are quite sensitive 
to the halo tail of the 6 He density used in the folding calculation and the IAS of 6 Li is likely to have 
a halo structure similar to that established for the ground state of 6 He. 

PACS numbers: 24.50.+g, 25.60.Bx, 25.60.Lg, 



The knowledge about the symmetry part of the nuclear 
equation of state (EOS) is vital for the understanding of 
the dynamics of supernovae explosion and the formation 
of neutron stars The symmetry part of the nuclear 

EOS is actually determined by the nuclear matter (NM) 
symmetry energy S(p) defined in terms of a Taylor series 
expansion of the NM binding energy B(p, S) as 



B(p,S) = B(p,0) + S(p)S 2 + O(5 4 ) 



(1) 



where 5 = (p n — p p )/p is the neutron-proton asymmetry 
parameter. The contribution of 0(<5 4 ) and higher-order 
terms in Eq. (|T|> . i.e., the deviation from the parabolic 
law was proven to be negligible 0, The NM sym- 
metry energy determined at the NM saturation density, 
E sym = S(po) with po « 0.17 fm -3 , is widely known 
in the literature as the symmetry energy or symmetry 
coefficient. Although numerous nuclear many-body cal- 
culations havepredicted E sym to be around 30 MeV (see, 
e.g., Refs. HH EH d), 

a direct experimental deter- 
mination of -E sym still remains a challenging task. One 
needs, therefore, to relate E syia to some experimentally 
inferrable quantity like the neutron skin in neutron-rich 
nuclei 0,0 IS El or the fragmentation data of heavy- ion 
(HI) collisions involving N £ Z nuclei [U [H El . An 
accurate estimate of the E syul value is also very impor- 
tant for the nuclear astrophysics. For example, a small 
variation of £ sym , used as input for the hydrodynamic 
simulation of supernovae, affects significantly the elec- 
tron capture rate in the "prompt" phase of type II su- 
pernovae [2|. Another example is a calculation of NM 
and masses of finite nuclei using Skyrme forces which 
shows that the neutron-rich NM does not collapse only if 
the corresponding E sym value is within the range 28 — 31 
MeV. -E S ym is also an important input for the study of the 
density dependence S(p) based on transport-model sim- 
ulation of the HI collisions (see Ref. [llj and references 



therein) , and the most recent transport-model results fa- 
vor E sym w 31 - 32 MeV [l3|. 

Within the frame of any microscopic model for asym- 
metric NM, the symmetry energy depends strongly on 
the isospin dependence of the nucleon-nucleon (NN) in- 
teraction used therein 0,0] ■ Therefore, the E sym value 
can be indirectly tested in a charge exchange (isospin- 
flip) reaction which has been known for decades as a good 
probe of the isospin dependence of the effective NN in- 
teraction ^4|. Although the isospin dependence of the 
optical potential (OP), known by now as Lane potential 
|l5j, has been studied since a long time, there has been a 
considerable interest recently in studying the isospin de- 
pendence of the OP in the quasi-elastic scattering reac- 
tions measured with unstable neutron-rich beams. Based 
on the isospin symmetry, the nucleus- nucleus OP can be 
written in terms of an isovector coupling 15] as 



U(R) = U (R) 



(2) 
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where t is the isospin of the projectile a and T is that of 
the target A. For a proton-induced scattering reaction, 
the second term of Eq. contributes to both the elastic 
(p, p) scattering and (p, n) charge exchange reaction |16| . 
While the relative contribution by the Lane potential U\ 
to the elastic (p,p) cross section is small and amounts 
only a few percent for a neutron-rich target |l7l Il8j , it 
determines entirely the (Fermi-type) AJ V = + tran- 
sition strength of the (p, n) reaction leading to an iso- 
baric analog state (IAS). Therefore, the (p,n) reaction 
has been so far the main tool in studying the isospin de- 
pendence of the proton- nucleus OP. Since this isospin de- 
pendence should be better tested in the charge exchange 
reactions induced by the neutron-rich beams, we consider 
in the present work thep( 6 He, 6 Li*)n reaction measured 
by Cortina-Gil et sd. [ljj with the secondary 6 He beam at 
Euh = 41. MeV. Given a large neutron-proton asym- 
metry (6 = 1/3) of the unstable 6 He nucleus, the mea- 
sured p( 6 He, 6 Li*)n cross section for the transition con- 
necting the ground state (g.s.) of 6 He (T = T z = 1) and 
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its isobaric analog partner (T = 1,T Z = 0, J v = + ex- 
cited state of 6 Li at 3.563 MeV) is indeed a good probe of 
the isovector coupling in the 6 He+p system. In the two- 
channel approximation, the elastic p( 6 He, 6 He)p scatter- 
ing and charge exchange p( 6 He, 6 Li*)n cross sections can 
be obtained from the solutions of the following coupled 
channels (CC) equations 0] 

[K p + U P (R) - E p ] Xp (R) = -^U^XniR), (3) 
[K n + U n {R) - E n ] X n{R) = --^tfi (4) 

Here Kpt n \ and Ep( n ) are the kinetic-energy operators 
and center-of-mass energies of the 6 He+p and 6 Li* + n 
channels, with the energy shift due to the Coulomb en- 
ergy and Q- value of the p( 6 He, 6 Li*)n reaction properly 
taken into account. U P (R) is the OP in the 6 He+p chan- 
nel and U n {R) is that in the 6 Li* + n channel. In addi- 
tion to the charge exchange jp( 6 He, 6 Li*)n and elas- 
tic p( 6 He, 6 He)p scattering |2y] data measured at 41. 6A 
MeV, a total reaction cross section <tr = 409 ±22 mb was 
also measured [2l| for the 6 He-t-p system at a slightly 
lower energy of 36A MeV. Thus, these data sets are the 
important constraints for the 6 He+p OP at the consid- 
ered energy. 

To link the Lane potential U\ to the isospin dependence 
of the NN interaction, we have used the folding model 
[rH fl8l ] to calculate Uq and U\ using the explicit pro- 
ton and neutron g.s. densities of 6 He and the CDM3Y6 
density- and isospin dependent NN interaction [22| ■ This 
interaction is based on the M3Y interaction fo(i)(s) de- 
duced from the G-matrix calculation 23] using the Paris 
NN potential, with the energy- and density dependences 
included explicitly as 

v 0(1) (E,p,s) = (1- 0.0026^0(1) (P)«0(i)(«). (5) 

where E is the bombarding energy (per projectile nu- 
cleon) , 



F o(i) (p) = Co(i) [1 + " exp(-/3p) - jp] , 



(6) 



and the explicit expression of the finite-range Wo(i)(s) in- 
teractions is given in Ref. |3( ■ Parameters of the isoscalar 
par t F (p) of the density dependence © were chosen 
|22| to reproduce saturation properties of the symmetric 
NM in the Hartree-Fock (HF) approximation and tested 
in the folding analyses [22, |24| of the elastic refractive 
nucleus-nucleus and a-nucleus scattering to infer realistic 
estimate for the nuclear incompressibility [K (po) « 250 
MeV]. In a similar manner, we try now to probe the 
isovector part F\ (p) in the CC analysis of the charge ex- 
change p( 6 He, 6 Li*)n reaction and, using the HF method 
of Ref. :3| , to accurately estimate the symmetry energy 
E sym . Since the CDM3Y6 interaction is real, the fold- 
ing model is used to calculate the real parts of the OP, 
Vo(i)=Re [Tom j which is further scaled by a complex fac- 



itly into account, one obtains from Eq. 

U P (R) = Vo(R)-^^ (Nn + iNj), (7) 
U n (R) = V (R)(N R + iN 1 ). (8) 

In the CC calculation, Z7 P („) are each added by a spin- 
orbital potential whose parameters were fixed as taken 
from the systematics by Becchetti and Greenless [25j . 
and U p is added further by a Coulomb potential between 
a point charge and a uniform charge distribution of ra- 
dius R c = 1.35 A 1 / 3 fm. The form factor (FF) of the 
p( 6 He, 6 Li*)n reaction, to be used in the right-hand sides 
of Eqs. © and (2), is given by Eq. J2J as 



Upn(R) = 



2y/N - Z 

A 



V2 



Ul(R) = yVi(fl)[l 



tor to obtained Uony By taking isospin coupling explic- 



i(Nj/N n )]. 

(9) 

Thus, the scaling factors JVr(i) of the real and imagi- 
nary parts of the OP are the main parameters to be 
determined from the CC description of the charge ex- 
change p( 6 He, 6 Li*)n and elastic scattering p( 6 He, 6 He)p 
data which should also be constrained by a total reaction 
cross section ctr ~ 400 mb (an empirical value expected 
at the considered energy |2j]). To have as less as possi- 
ble free parameters, we have used a simple assumption 
(JSJ to scale the real folded (p, n) FF by the same relative 
complex strength as that used in the elastic channel [26| . 
Our only nuclear structure input is the 6 He g . s . density 
and we have considered in this work two different choices: 
the microscopic density given by the cluster-orbital shell 
model approximation (COSMA) |27l |28| and that ob- 
tained recently [2!| based on the independent particle 
model (IPM). The CC calculation was done with the non- 
relativistic code FRESCO j^D] using the inputs for mass 
numbers and incident energies given by the relativisti- 
cally corrected kinematics |3l|. For a checking purpose, 
the CC results plotted in Fig. ^ were also compared with 
those given by the code ECIS97 1^2 (which takes exactly 
into account the relativistic kinematics) and the two sets 
of calculated p( 6 He, 6 Li*)n cross sections turned out to 
be nearly identical. 

We found that Nr « 0.85 and iVj w 0.55 which were 
mainly determined by the fit to the elastic scattering data 
and by a constraint that the calculated cjr. is around 400 
mb. It is noteworthy that a continuum-discretized cou- 
pled channels (CDCC) calculation of the elastic 6 Hc+p 
scattering at about the same energy by Mackintosh and 
Rusek |3J3] has shown that the real dynamic polariza- 
tion potential due to 6 He breakup is repulsive in the cen- 
ter and at the surface, so that a renormalization factor 
2Vr < 1 of the real folded OP is well expected. The rel- 
ative strength Ni/Nr « 0.65 (used further to calculate 
FF for the p( 6 He, 6 Li*)n reaction) also agrees reason- 
ably with the CDCC results |3^] which give the ratio of 
volume integrals of the imaginary and real parts of the 
6 He+p OP around 0.6. 

The p( 6 He, 6 Li*)n cross section given by the COSMA 
density [28j was found to have a shape very close to that 
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FIG. 1: CC results for the charge exchange p( 6 He, 6 Li*)n 
cross section at £i a b = 4 1.6 A MeV in comparison with the 
data measured by Cortina-Gil et al. |Til . 



of the measured angular distribution (see Fig.^). Since 
the strength Ni/Nr was fixed at 0.65, the CC description 
of the p( 6 He, 6 Li*)n data could be improved only by fine 
tuning the strength C\ of the isovector part of the den- 
sity dependence © of the CDM3Y6 interaction. One 
can see that the best fit is achieved when C\ is about 
10% stronger than the isoscalar strength Co. We have 
further performed the HF calculation of asymmetric NM 
with the same isospin- and density dependent CDM3Y6 
interaction using the method described in Ref. 3]. The 
density dependence of the NM symmetry energy S(p) 
obtained with the same isovector strengths C\ as those 
used in Fig. ^ is shown in Fig. |2J and one can deduce 
31 ± 1 MeV from our HF results. This 



sym 



easily E, 

result should be complementary to the nuclear structure 
studies which relate the slope of the EOS of asymmetric 
NM and the associated E sym value to the neutron skin, a 
method first suggested by Alex Brown 0. If one adopts, 
e.g., a neutron-skin AR w 0.1 — 0.2 fm for 208 Pb then 
a systematics based on the mean-field calculations (see 
Fig. 7 of Ref. 0) gives E sym « 27 - 31 MeV (this value 
is compared with our HF result in Fig. 0). The main 
methods to determine the neutron skin are either the 
analyses of elastic (p, p) scattering on stable N ^ Z tar- 
gets 0, |35| or studies of asymmetric NM and structure 
of finite nuclei [Ioll36| . However, the uncertainty remains 
still rather high and AR for 208 Pb nucleus was found to 
be ranging from 0.083-0.11 fm [H| to 0.13 ± 0.03 fm [3 
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FIG. 2: Density dependence of the NM symmetry energy S(p) 
predicted by the HF formalism of Ref. |3] using the same 
isovector strengths Ci as those used in Fig. in comparison 
with the empirical points deduced from the neutron-skin 
and HI fragmentation |T^| data. 



or around 0.17 fm |3J, and up to about 0.22 fm [36j. A 
more accurate determination of AR is expected from the 
measurement of parity-violating electron scattering 
and it might be used for a more precise determination of 
E S ym- Our result is also complementary to the nuclear 
reaction studies based on the transport-model simula- 
tions fill IT^ . For example, a recent antisymmetrized 
molecular dynamics (AMD) analysis of the HI fragmen- 
tation data [12 obtained S(p « 0.08 fm~ 3 ) w 18-20 
MeV (see Fig. 01 which gives a useful information on 
the low-density part of S(p). Since the charge exchange 
p( 6 He, 6 Li*)n cross section is peaked at the most for- 
ward angles (see Fig. 0), our CC analysis has probed 
mainly the surface part of the FF which, in turn, is de- 
termined by the low-density part of v\ (E, p, s). The fact 
that our HF calculation reproduces quite well the empir- 
ical half-density point of S(p) ^3 shows the reliability 
of the isospin dependence of the CDM3Y6 interaction. 
Note that the Gogny-AS effective NN interaction used 

31 



J sym 



in the AMD calculation of Ref. [12j also gives E. 
MeV at po (see Fig. 1 of Ref. |23), in a very close agree- 
ment with our HF result. Thus, the HF result shown in 
Fig. [3 should provide a realistic description of the EOS 
of asymmetric NM for densities up to about 2 po. 

Besides the conclusion on the symmetry energy, we 
have found further that the measured p( 6 He, 6 Li*)n data 
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FIG. 3: CC results for the elastic scattering p( 6 He, 6 He)p and 
charge exchange p( 6 He, 6 Li*)n cross sections at i?i a b = 41.6^4 
MeV, given by the OP and FF obtained with 3 choices of 
the 6 He g . s . density (see text), in c omp arison with the data 
measured by Cortina-Gil et al. |19ll20|. 
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FIG. 4: DWBA results for the inelastic p( 6 He, 6 He*)p' scat- 
tering to the 1.87 MeV 2+ state of 6 He at _Ei ab = 40.9A MeV, 
given by the OP and FF obtained with 3 choices of the 6 He g . s . 
density (see text), in comparison with the data measured by 
Lagoyannis et al. |4l|. 



are also sensitive to the halo tail of 6 He nucleus. Since the 
IAS of 6 Li is just an isobaric analog partner of 6 He g s , the 
isospin symmetry implies |30l | that this IAS of 6 Li should 
have about the same halo structure as that of 6 He g . s ., and 
the central OP's in the entrance and exit channels should 
be well described by Eqs. 0-©- To study this effect, 
we have used in the folding calculation two versions of 
the COSMA density jH which have RMS radii of 2.57 
and 2.48 fm as well as the IPM density [2£j which has the 
RMS radius of 2.46 fm. Although the IPM density was 
shown [2^ to give a good description of the interaction 
cross section measured with 6 He beam at high energies, it 
was calculated in the independent particle model which 
does not account for the di-neutron correlation in 6 He. 
In this sense, the COSMA densities are more accurate 
and should give a better description of the p( 6 He, 6 Li*)n 
data if the IAS of 6 Li has the same halo structure as 
6 He g . s .. The CC results obtained with 3 different choices 
of the 6 He g . s . density are shown in Fig. [21 and one can 
see that the COSMA density is indeed more appropriate 
than IPM density and gives consistently good description 
of both the elastic scattering and charge exchange data. 
With the same iVR(i) factors used throughout these cal- 
culations, the COSMA densities with RMS = 2.57 and 
2.48 fm give ctr = 408 and 399 mb, respectively, quite 
close to the empirical value of about 400 mb. The IPM 



density gives or, = 382 mb which is slightly smaller than 
400 mb. The description of the elastic scattering and 
charge exchange data by the IPM density can only be 
improved by further reducing Nj factor but then <tr be- 
comes significantly smaller than the empirical value. One 
can see in Fig. [3] that from two versions of the COSMA 
density the p( 6 He, 6 Li*)n data seem to favor that with 
RMS=2.57 fm. Such a RMS radius agrees with that pre- 
dicted earlier by the three-body calculation of 6 He which 
takes into account the dynamic correlation between the 
a-core and di-neutron |40| . 

In addition to the elastic p( 6 He, 6 He)p scattering 
and charge exchange p( 6 He, 6 Li*)n reaction, inelastic 
p( 6 He, 6 He*y scattering to the 1.87 MeV 2+ state of 
6 He has also been measured 0] at a nearby energy of 
40. 9A MeV. This (unbound) excitation of 6 He has been 
a subject of extensive analyses of the elastic and inelas- 
tic 6 He+p scattering either in the distorted wave Born 
approximation (DWBA) or CC formalism. In particular, 
the halo effect has been shown to be significant in the 
2 + inelastic scattering channel (see Fig. 3 of Ref. |41|). 
Therefore, a folding analysis of the 2 + inelastic scattering 
data using the same 6 He+p OP's as those used in Fig- El 
would be quite complementary to the results discussed 
above for the charge exchange reaction. The real 2 + in- 
elastic FF was calculated by the folding model 0, 0] 
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using a simple ansatz for the transition density, where the 
proton and neutron parts of the ( 6 He g s . — » 6 He* + ) tran- 
sition density are given by deforming proton and neu- 
tron parts of the 6 He g . s . density with the deformation 

lengths 8$ an d 8%+ determined recently in a CC anal- 
ysis of the 4 He( 6 He, 6 He) 4 He reaction The folded 
2 + inelastic FF was further scaled by the same complex 
factor (1 + iNi/Nr) as that used to obtain the charge 
exchange FF in Eq. ©. The DWBA results obtained 
with three choices of the 6 He g . s . density are plotted in 
Fig. 0] and one can see that the best description of the 
p( 6 He, 6 He*)p' data is again given by the COSMA den- 
sity. The deficiency of the IPM density is about the same 
as that shown above in the calculated p( 6 He, 6 Li*)n cross 
sections. 

Finally, we note that a folding analysis of the present 
charge exchange p( 6 He, 6 Li*)n data has been done ear- 
lier |19l l2l| using the famous JLM (complex) G-matrix 
interaction. Although the real, imaginary and isovector 



strengths of the JLM interaction were adjusted to the 
best fit of the elastic scattering and charge exchange data, 
these analyses seem to be unable to give a good descrip- 
tion of the last data points of the measured p( 6 He, 6 Li*)n 
cross section no matter what density distribution of 6 He 
is used in the folding calculation (see, e.g., Fig. 2 of 
Ref. [HI)- Our CC results represent, therefore, an ac- 
curate alternative description which also provides impor- 
tant input for the description of the EOS of asymmetric 
NM. The future measurements of the charge exchange 
reactions induced by the neutron-rich beams would be 
very valuable in studying the isospin dependence of the 
nucleus-nucleus interaction and making a more reliable 
conclusion on the symmetry energy. 
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